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Section  1 
INTRODUCTION 


The  axial-flow  fans  and  compressors  used  in  aircraft  gas-turbine 
propulsion  systems  frequently  operate  in  a  range  where  the  flow  over  the  blades 
is  transonic.  This  situation  arises  at  high  values  of  the  rotor  angular  ve¬ 
locity,  where  the  resultant  of  the  axial  inlet  speed  of  the  gas  and  the 
rotational  speed  of  the  blades  reaches  the  speed  of  sound  at  some  station  along 
the  blade  span.  Under  these  circumstances,  the  moving  blade  experiences  an 
inlet  relative  flow  that  is  supersonic  over  the  outer  portion  of  the  blade, 
and  is  subsonic  near  the  hub.  Shock  waves  are  formed  at  the  leading  edge  of 
the  outer  portion;  these  propagate  upstream  through  the  subsonic  axial  inlet 
flow,  and  are  radiated  away  from  the  inlet  itself. 

The  objective  of  the  research  reported  here  was  to  develop  a  method 
for  incorporating  these  radiating  waves  in  the  numerical  codes  that  are  used 
for  calculations  of  such  flows.  Calculations  of  this  type  have  come  into 
increasingly  heavy  use  in  recent  years  as  a  design  tool,  for  screening  candidate 
configurations.  However,  the  calculation  methods  used  at  present  ignore  the 
effects  due  to  the  radiating  waves,  and  assume  instead  that  the  flow  at  the 
inlet  is  uniform,  and  adequately  described  by  the  two-dimensional  flow  approxi¬ 
mations  at  any  radial  location. 

The  adequacy  of  these  approximations  has  been  the  subject  of  con¬ 
siderable  study  in  recent  years,  especially  because  of  the  transonic  nature 

of  the  flow.  If,  for  example,  the  wave  pattern  radiated  by  the  blades  is 

1  2 

estimated  by  the  two-dimensional  simple-wave  model,  ’  then  the  flow  has 


T!  Rae,  W.J.  and  Lordi,  J.A.,  "A  Study  of  Inlet  Conditions  for  Three- 

Dimensional  Transonic  Compressor  Flows",  Calspan  Report  No.  XE-6129-A-4 
(June  1978). 

2.  Fink,  M. R.  "Shock  Wave  Behavior  in  Transonic  Compressor  Noise  Generation" 
Trans  ASME  93  (D)  Journal  of  Engineering  for  Power  No.  4  (October 

1971)  397-403. 
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singular  behavior  at  the  sonic  circle,  i.e.,  at  the  radial  station  where  the 

inlet  relative  Mach  number  is  equal  to  one.  Moreover,  the  distribution  of 

inlet  mass  flow  over  the  supersonic-inlet  portion  is  conventionally  taken  as 

3 

that  given  by  the  unique- incidence  condition,  applied  at  each  radius.  This 
approximation  neglects  the  radial  interactions  in  the  flow,  and  the  resulting 
redistribution  of  mass  flow  that  is  known  to  occur. 

In  addition  to  the  inadequacy  of  these  approximations,  it  is  reasonable 
to  expect  that  the  absence  of  a  radiation  condition  would  introduce  numerical 
anomalies  at  the  grid  boundaries.  The  types  of  anomalous  wave  patterns  ex¬ 
pected  are  shown  in  Figure  1,  which  is  a  section  through  the  blade  row  at  a 
constant  radius,  outside  the  sonic  circle.  This  figure,  taken  from  a  previous 
NAVAIR- sponsored  study  at  Calspan,*  suggests  that  waves  which  propagate  upstream 
from  the  blade  row  would  be  reflected  from  the  upstream  grid  boundary.  As 
pointed  out  in  Reference  l,  previous  attempts  to  calculate  a  supersonic-tip- 
speed  case  had  encountered  wavelike  motions  at  the  radial  stations  beyond  the 
sonic  circle.  These  motions  were  assumed  to  be  caused  by  the  lack  of  a  radi¬ 
ation  condition,  and  the  objective  of  Reference  1,  as  well  as  the  present  work, 
was  to  apply  a  radiation  condition  in  order  to  allow  the  waves  to  escape,  and 
to  learn  from  the  resulting  solutions  how  to  simplify  the  procedure. 

The  approach  that  was  followed  in  Reference  1  was  to  use  the  linear- 

theory  formulation  of  the  radiation  condition,  in  conjunction  with  a  computer 

code  which  solves  the  three-dimensional  transonic  flow  in  the  nonlinear  small- 

4  5 

disturbance  approximation.  ’  The  two  solutions  were  matched  along  a  plane 


3.  Lichtfuss,  H.J.  and  Starken,  H.  "Supersonic  Cascade  Flow",  pp.  37-149 
of  Progress  in  Aerospace  Sciences  Vol.  15  ed.  by  D.  Kuchemann  Pergamon 
Press,  New  York  (1974). 

4.  Rae,  W.J.,  "Computer  Program  for  Relaxation  Solutions  of  the  Non-Linear 
Small-Disturbance  Equations  for  Transonic  Flow  in  an  Axial  Compressor 
Blade  Row",  AFOSR-TR-78-0855,  (April  1978). 

5.  Rae,  W.J.  "Calculations  of  Three-Dimensional  Transonic  Compressor  Flow 
Fields  By  A  Relaxation  Method"  Journal  of  Energy  1  pp.  284-296 
(1977). 
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Figure  1  VJRvB  PRTTERMS  EXPECTED 


located  a  chord  length  or  so  upstream  of  the  blade  row;  the  numerical  solution, 
at  a  given  stage  in  the  iterations,  was  used  to  calculate  the  amplitude  coef¬ 
ficients  in  the  Fourier- Bessel  series  representation  of  the  outgoing  waves 
in  the  linear  theory-  The  velocity  field  derived  from  this  series  was  then 
used  as  the  boundary  value  for  the  next  step  of  the  numerical  solution,  and 
the  process  was  repeated  at  regular  intervals. 

In  Reference  1,  the  case  chosen  for  application  of  this  technique 
was  a  heavily  loaded  rotor  with  a  relative  Mach  number  of  1.3  at  the  tip.  A 
series  of  attempts  to  calculate  this  flow,  using  the  uniform-inlet  boundary 
conditions,  had  led  to  solutions  in  which  wavelike  disturbances  propagated 
back  and  forth  through  the  flow.  It  was  expected  that  introduction  of  the 
radiation  condition  would  remove  these  wavelike  oscillations.  However, 
it  did  not;  many  attempts  (reported  in  Reference  1)  were  made,  with  adjustments 
in  the  frequency  of  applying  the  radiation  condition,  but  none  of  these  succeeded 
in  removing  the  waves.  The  work  of  Reference  1  concluded  with  the  speculation 
that  the  case  attempted  was  too  difficult  a  test  for  the  method,  in  that  the 
blades  were  quite  heavily  loaded,  and  the  initial  flowfield  too  filled  with 
steep  waves. 

Accordingly,  the  present  investigation  was  undertaken,  in  which  a 
new  calculation  was  started  from  the  very  beginning.  The  blade  geometry  and 
operating  conditions  were  chosen  so  as  to  have  relatively  modest  values  of  the 
loading  and  Mach  number.  For  this  case,  the  procedure  described  above  did  lead 
to  a  converged  solution. 

However,  the  calculation  was  then  repeated  with  the  uniform- inlet 
boundary  condition,  in  order  to  verify  that  the  radiation  condition  was  the 
factor  responsible  for  the  convergence.  It  was  found  that  this  calculation 
also  converged,  and  converged  to  essentially  the  same  answer  as  that  found 
using  the  radiation  condition. 

To  explore  the  problem  further,  the  heavily  loaded  case  was  then 
calculated  again  -  once  with  the  radiation  condition  and  once  with  the 
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uniform- inlet  condition.  Both  of  these  calculations  were  started  from  an 
initial  condition  of  uniform  flow,  and  both  led  to  the  same  phenomenon  ob¬ 
served  in  the  earlier  work,  namely,  a  series  of  waves  was  found  to  oscillate 
back  and  fortr  in  the  flow  field.  The  detailed  distributions  of  these 
oscillations  at  any  stage  of  the  iteration  process  were  virtually  the  same 
for  either  of  the  boundary  conditions  used  at  the  inlet. 

Thus  the  final  conclusion  reached  is  a  negative  one:  for  the  cal¬ 
culation  that  converged,  the  use  of  a  radiation  condition  led  to  results  which 
differed  negligibly  from  those  found  using  the  uniform- inlet  condition.  Toward 
the  end  of  the  work  reported  here,  a  similar  conclusion  was  reached  in  a  study 
done  at  the  MIT  Gas  Turbine  Laboratory  (Reference  6),  using  essentially  the 
same  formulation  of  the  radiation  condition,  in  conjunction  with  a  numerical 
code  that  solves  the  full  Euler  equations. 

The  steps  followed  in  reaching  this  conclusion  are  outlined  in  the 
Sections  below:  Section  2  contains  a  brief  review  of  the  numerical  method, 
including  some  improvements  that  were  made  in  the  basic  non-linear  small-dis¬ 
turbance  program  while  Section  3  describes  the  present  calculations.  Section 
4  contains  a  discussion  of  the  conclusions  reached. 

This  report  relies  heavily  on  Reference  1,  and  is  intended  to  be 
read  in  conjunction  with  it.  Much  of  the  notation,  blade-geometry  description, 
and  computer-program  details  mentioned  here  are  taken  directly  from  Reference 
1;  those  details  have  not  been  repeated  in  this  text. 


6.  Berenfeld,  S. ,  "A  Study  of  Boundary  Conditions  for  Transonic  Compressor 
Flows",  M  Sc.  Thesis,  Massachusetts  Institute  of  Technology 
(June  1979). 


After  Reference  4  was  published,  it  was  found  that  an  extended  form 
of  the  Thomas  algorithm  could  be  used,  in  which  the  periodic  boundary  condi¬ 
tions  are  enforced  at  every  stage  of  the  relaxation  process,  rather  than 
being  approached  only  as  the  solution  converges.  This  extended  form  is 
described  in  Reference  9. 


The  details  of  the  method  are  as  follows:  the  grid  points  along 
the  line  being  relaxed  are  numbered  L  =  1,  2,  3,  ...,  LMX;  the  difference 
equation  has  the  form  of  a  three-term  recurrence  relation: 


+  BL4>L  +  cV  '  =  DL  L  =  2,  3,  ...,  LMX 


where  $  is  the  velocity  potential,  and  the  coefficients  RL ,  ...»  DL  are  known 

quantities  at  each  value  of  L.  The  periodicity  condition  to  be  enforced  up¬ 
stream  and  downstream  of  the  blades  is 


.  /  .LMX 

<j>  =  q>  +  G  A  (p 


where  6*0  upstream  of  the  blades,  G  *  1  downstream  of  the  blades,  and  A 
is  the  potential  jump  at  the  trailing  edge. 


The  extended  Thomas  algorithm  is  implemented  as  follows:  let 


£  (»>  -  O  ,  F(1)  =  o  ,  So)  *  1 


Then  calculate 


4  4  /  [  C  ‘  £  (  L  -  t  )  *  B  L  ] 


D‘  CL  r  (  l  -  f)  3  /  1C  E(L-r)  *  B  L  -  2,3,...,  /.MX 


IC  £  <’*.  0*8  J 


»'  'v**i  si f  x  .  and  !ial«h.  .J.l  .  ,  The  Theory  of  Splines 

, ...  -  'j-j  .  iti'rt.  »,  n  Ve*  York  ( 1 9f< . )  pp .  14-16. 


I  •#«  •  #  #•  •  •  •  •  •  ■ 


Next ,  set 


T  (LM  X)  -  I  ,  V^MX)  ■=  0 


and  calculate 

T(L)  -  E(L)  T  (Lt- 1)  *■  S(L) 

V(L)  -  E  CL)  V  (L  +  t)  *  FLL  ) 
Then  the  potential  at  LMX  is  given  by 


L 


LMX  -  }  ,  LMX  -  2> 


I- MX. 


ftL”*LV(Z)  -  £A<f)  ]  -  CLnXVUMX  -i)  +  DLH* 


R*'”*  TCZ)  +  T  (LMX-  1)  +  BLMX 


Z 


Once  this  is  known,  the  other  values  of  LMX  are  found  from  the  recursion  re¬ 
lations: 


<t>L  =  E(L)  (f)L*1  +  F(L)  +  SCL)  (pLr1*,  L  -  2,3,...,  LMx-i 
<p‘  -  <t>LM*  +  £A<p 

At  points  on  the  blade,  the  normal  Thomas  algorithm  is  used,  as 
before,  by  setting  S,  T,  and  V  equal  to  zero. 

Correction  Form 


The  accuracy  of  the  solution  can  be  improved  by  using  the  "correction" 
form  of  the  algorithm,  i.e.,  on  each  iteration  the  change  in  the  velocity  po¬ 
tential  is  calculated: 

C  *  <f>  -  <f> 

+ 

rather  than  the  new  value  of  the  potential  itself,  <j>  .  This  form  minimizes 
the  arithmetic  problems  that  might  arise  from  taking  small  differences  of  large 
numbers . 
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9  •  ■>  #«< 
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Exit-Plane  Mach  Number 


The  final  change  made  in  the  relaxation  method  deals  with  the  Mach 
number  variation  at  the  downstream  edge  of  the  grid.  At  that  location,  the 
velocity  component  is  constrained  to  equal  a  certain  constant,  dictated  by 
the  requirement  of  mass-flow  conservation.  In  the  previous  version  of  the 
program,  this  requirement  was  met  by  setting  the  potential  values  at  the  down¬ 
stream  edge  of  the  grid  equal  to  the  value  corresponding  to  a  parabola  through 
the  two  stations  just  upstream,  and  having  the  required  slope.  The  revision  in 
corporated  in  the  present  work  is  to  set  the  velocity  equal  to  the  desired 
value  at  the  next-to-last  station,  and  to  use  this  velocity  in  solving  for 
the  potential  there.  The  potential  at  the  last  station  is  then  found  by  a 
linear  extrapolation,  using  the  desired  slope.  A  similar  procedure  was  fol¬ 
lowed  by  Collins  and  Krupp.10 


10.  Collins,  D.J.,  and  Krupp,  J.A.,  "Experimental  and  Theoretical  Investiga 
tions  in  Two-Dimensional  Transonic  Flow",  AIAA  Journal  12^  (1974) 
pp.  771-778. 


Section  3 

RESULTS  OF  PRESENT  CALCULATIONS 


Blade-Row  Geometry 

As  noted  in  the  Introduction,  the  case  chosen  for  the  previous  work 
was  felt  to  be  too  stringent  a  test  of  the  radiation  condition,  because  of 
the  high  degree  of  loading  and  because  the  initial  solution  that  was  used 
already  contained  steep  wavelike  disturbances.  Accordingly,  the  blade  geometry 
and  operating  conditions  for  the  present  calculations  were  carefully  chosen 
so  as  to  avoid  these  problems.  Flat-plate  blades  were  used,  in  order  to 
eliminate  blade  thickness  as  a  possible  source  of  choking,  and  the  blades 
were  twisted  so  as  to  present  a  constant  angle  of  incidence  of  six  degrees  to 
the  relative  inlet  flow.  This  incidence  was  sufficient  to  insure  subsonic 
flow  over  most  of  the  region  downstream  of  the  blades,  thus  minimizing  wave 
reflection  at  the  downstream  edge  of  the  grid  as  a  mechanism  for  instability 
in  the  calcul  '.tions.  The  blade  row  had  30  blades,  with  a  hub/tip  radius  ratio 
of  0.75  and  a  solidity  C*  /  LT  equal  to  0.5.  The  axial  Mach  number  was  set 
equal  to  0.5,  and  the  Totational  Mach  number  at  the  tip  was  0.980,  giving  a 
distribution  of  inlet  relative  Mach  number  that  varied  from  0.889  at  the  hub 
to  1.100  at  the  tip. 

The  grid  had  60  points  in  the  axial  direction,  20  in  the  blade-to- 
blade  direction,  and  10  radially.  The  grid  extended  from  two  axial  chords 
upstream  of  the  leading  edge  to  two  axial  chords  downstream  of  the  trailing 
edge,  and  was  stretched,  as  described  in  Reference  4,  to  concentrate  points 
near  midchord.  This  grid  had  the  optimum  grid  size  ratio  for  shock-wave 
capturing  at  mid  span  (see  Reference  4) . 

Iteration  Sequence 

The  Fourier- Bessel  series  was  calculated  using  11  terms  in  the 
***  -series  (•**■  0  through  10)  and  10  terms  in  the  A  -series  ( «  1  through 
10).  On  the  first  run,  the  quantities  ^m*|,  and  N^b^ (see  Section  4 


of  Reference  1)  were  calculated  and  stored  on  tape.  These  values  are  shown 
in  Table  1,  where  the  dotted  line  separates  the  inodes  that  propagate  (above 
the  line)  from  those  that  are  damped  (below  the  line).  The  calculations  then 
proceeded  through  a  series  of  sixteen  radial  iterations,  each  of  which  used 
60  axial  iterations.  Thus  the  solution  at  each  grid  point  was  iterated  960 
times.  The  fifth  axial  grid  station,  lying  at  0.9752  axial  chordlengths  up¬ 
stream  of  the  leading-edge  plane,  was  chosen  as  the  matching  plane.  For  the 
first  60  iterations  (i.e.,  for  the  first  radial  sweep)  uniform-flow  conditions 
were  imposed  at  the  inlet,  and  thereafter  the  Fourier- Bessel  series  was  used, 
with  the  amplitude  coefficients  re-evaluated  every  60  iterations  (i.e., 

prior  to  every  radial  iteration).  The  angle  of  incidence  and  the  relaxation 
factors  used  were  increased  gradually:  incidence  of  2°  for  the  first  240 
iterations,  4°  for  the  next  60,  and  6°  thereafter,  while  the  relaxation 
factors  (for  supersonic  and  subsonic  points)  were  both  taken  as  0.2  for  the 
first  120  iterations,  0.5  for  the  next  240,  and  0.8  thereafter.  An  attempt 
to  raise  the  subsonic  relaxation  factor  to  1.2,  starting  at  the  601st  itera¬ 
tion,  led  to  a  divergent  solution. 

The  convergence  of  the  solution  can  be  seen  in  the  history  of  the 

local  Mach  number  at  a  fixed  point.  Figure  3  shows  this  quantity  at  t,  *  O  , 

X/  =  0.0126,  p  -  1.7419.  This  point  is  located  on  the  suction  side, 

just  downstream  of  the  leading  edge,  and  at  a  radius  where  the  inlet  relative 

Mach  number  is  1.0043.  Thus  it  is  in  a  region  of  peak  acceleration.  It  is 
clear  that  the  solution  has  essentially  converged,  at  this  stage  of  the  itera¬ 
tions.  Another  measure  of  the  convergence  is  the  average  residual,  which 
after  960  iterations  had  decreased  to  4.2%  of  the  sum  of  the  absolute  values 
of  the  five  terms  in  the  finite-difference  equation.  Previous  experience 
with  this  program  indicates  that  further  reductions  of  the  residuals  would 
require  many  more  iterations,  but  would  yield  insignificant  changes  in  the 
velocity  distributions. 

The  amplitude  coefficients  in  the  Fourier-Bessel  series  did  not  show 
the  same  degree  of  convergence.  Some  of  these  coefficients  were  changing  only 
slightly  at  the  end  of  the  calculation,  but  others  were  still  growing;  in 
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particular,  the  constant  term  f)go  was  still  increasing.  Table  2  shows  the 
values  of  these  coefficients  used  for  the  final  60  iterations.  It  is  obvious 
that  only  the  first  few  modes  are  present  in  the  numerical  data. 

Sonic- line  contours  from  the  final  iteration  are  shown  in  Figure  4; 
the  significance  of  these  results  is  discussed  at  greater  length  in  Section  4. 

Check  Runs 


After  the  solution  described  above  was  completed,  several  shorter 
runs  were  made,  in  order  to  verify  that  the  radiation  condition  was  in  fact 
the  critical  element  that  had  allowed  convergence.  In  these  shorter  runs, 
the  final  120  iterations  were  repeated,  without  the  radiation  condition,  and 
using  several  values  of  the  relaxation  factors.  These  runs  produced  the  sur¬ 
prising  result  that  the  solutions  were  virtually  identical  to  that  described 
above.  It  is  possible  that  120  iterations  is  insufficient  to  allow  instabil¬ 
ities  in  the  flow  to  develop,  but  it  seems  unlikely,  especially  in  view  of 
the  high  degree  of  identity  between  solutions  with  and  without  the  radiation 
conditions:  there  were  no  differences  between  the  two  solutions,  to  three- 
figure  accuracy.  The  fact  that  the  detailed  history  of  these  two  solutions 
was  so  nearly  identical  prompted  several  checks  of  the  program  logic,  to  verify 
that  the  radiation  condition  was  actually  being  used. 

At  this  point,  it  was  decided  to  repeat  the  calculation  of  the  highly 
loaded  case  done  in  the  previous  study.  This  was  done  to  find  out  whether  the 
changes  that  had  been  made  to  the  program  might  have  removed  some  previously 
undetected  errors  that  were  causing  the  oscillations  in  the  solution. 

Thus,  the  23-blade  heavily-loaded-rotor  case  described  in  the  previous 
report  was  repeated,  starting  from  the  very  beginning,  i.e. ,  from  a  uniform- 
flow  condition  throughout  the  flow  field.  Relaxation  factors  of  0.8  were  used 
at  all  points,  subsonic  and  supersonic,  throughout  the  entire  series  of  itera¬ 
tions.  To  start  the  calculation,  two  radial  sweeps  were  done,  each  with 
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60  axial  sweeps,  for  a  total  of  120  updates  of  the  solution  at  each  grid  point. 
Thereafter,  this  initial  run  was  continued  on  two  parallel  paths  -  with  the 
radiation  condition  and  with  the  uniform-flow  inlet  condition.  Both  of  these 
calculations  were  run  for  10  further  radial  iterations,  at  60  axial  iterations 
each,  bringing  the  total  number  of  iterations  to  720.  The  radiation  condition 
was  updated  every  60  iterations. 

Both  of  these  solutions  showed  the  same  type  of  oscillatory  behavior 
as  that  observed  during  the  previous  work.*  Moreover,  the  oscillating  time 
histories  of  the  potential  distribution  calculated  with  the  two  inlet  condi¬ 
tions  were  virtually  identical  to  each  other,  iteration  by  iteration. 


Section  4 
CONCLUSIONS 


The  conclusion  that  emerges  from  these  calculations  is  that  the 
boundary  conditions  used  at  the  inlet  have  a  negligible  effect  on  the  flow 
field  solution,  at  least  for  the  set  of  conditions  represented  by  the  case 
which  converged.  This  conclusion  is  in  accord  with  the  results  of  an  inde¬ 
pendent  study  recently  completed  at  the  MIT  Gas  Turbine  Laboratory.6  While 
differing  in  details,  the  essential  conclusion  reached  in  that  work  is  the 
same  as  that  reached  here. 

The  limits  over  which  this  conclusion  applies  are  not  established 
by  the  results  of  the  present  study,  since  only  a  single  case  has  been  run 
to  convergence.  For  the  case  which  did  not  converge,  it  is  true  that  the  inlet 
condition  had  no  effect  on  the  solution,  but  this  leaves  open  the  question 
whether  such  an  effect  was  suppressed  by  the  presence  of  the  oscillating  waves. 

The  original  hypothesis  on  which  this  work  was  based  was  that  the 
oscillations  in  the  solution  were  caused  by  wave  reflections  from  the  upstream 
boundary.  These  reflections,  in  turn,  were  assumed  to  be  generated  by  the 
imposition  of  the  uniform- inlet  conditions.  The  converged  results  from  the 
twisted  flat-plate  case  make  it  clear  that  the  latter  assumption  is  not 
always  realized;  the  Mach-number  contours  in  Figure  4  are  not  affected  by  the 
imposition  of  uniform-flow  conditions  at  the  matching  plane.  Thus  the  mech¬ 
anism  for  generating  oscillations  in  the  solution  is  more  complex;  it  may  be 
dependent  on  the  amplitude  of  the  disturbances  reaching  the  matching  plane, 
on  the  relaxation  factors  being  used,  and  on  other  sources,  that  are  unknown 
at  present. 


Several  aspects  of  the  calculations  reported  here  require  further 
elaboration:  the  number  of  iterations  used,  the  occurrence  of  supersonic  out¬ 
flow,  and  the  schedule  used  for  the  relaxation  factors. 


Number  of  Iterations 


The  magnitude  of  the  residuals  remaining  after  960  iterations  (for 
the  twisted  flat-plate  case  described  in  Section  3)  is  somewhat  larger  than 
the  standard  usually  considered  acceptable  for  a  finite-difference  solution. 
However,  it  is  also  characteristic  of  line-relaxation  methods  that  the  short- 
wavelength  errors  decay  very  slowly.  And  while  these  errors  may  yield  large 
values  of  the  residual,  they  do  not  have  a  major  effect  on  the  gross  flow  field, 

nor  do  they  lead  to  divergent  behavior  as  the  solution  progresses.  The  general 
nature  of  line-relaxation  methods  indicates  that  these  solutions  were  carried 
far  enough  to  provide  the  required  opportunity  for  the  radiation  condition  to 
make  its  presence  felt. 

Supersonic  Outflow 

The  angle  of  incidence  used  for  the  twisted  flat-plate  case  was 
raised  to  six  degrees  in  order  to  insure  that  the  flow  at  the  exit  plane  would 
be  subsonic.  As  noted  earlier,  this  was  done  to  remove  wave  reflections 
from  the  downstream  edge  of  the  grid  as  a  possible  source  of  oscillatory 
behavior.  The  Mach  number  distributions  reached  on  the  960th  iteration 
fell  just  short  of  this  goal,  in  that  the  exit-plane  Mach  number  at  the  radial 
station  next  to  the  outer  wall  was  1.003.  However,  there  was  no  evidence  of 
wave  reflections  in  this  vicinity. 

Relaxation-Factor  Schedule 


One  noteworthy  difference  between  the  twisted  flat-plate  case  and  the 
heavily- loaded-rotor  case  is  the  schedule  of  relaxation  factors  used.  These 
were  varied  from  0.2  to  0.8  in  the  former  case,  but  were  set  equal  to  0.8 
throughout  the  latter  case.  This  raises  the  question  whether  the  latter  case 
might  have  converged,  if  the  relaxation  factors  had  been  increased  gradually. 


Finally,  it  is  appropriate  to  return  to  consideration  of  the  longer- 
range  goal  of  this  research,  namely  to  shed  light  on  the  nature  of  the  subsonic/ 
supersonic  transition  that  occurs  near  the  sonic  circle.  The  twisted  flat-plate 
solution  does  not  offer  a  great  deal  of  information  concerning  this  question. 
Specifically,  the  sonic-line  contours  found  in  the  vicinity  of  the  sonic  circle 
do  not  suggest  any  simple  interpretation;  they  are  not  aligned  at  any  recog¬ 
nizable  Mach  angle,  and  the  radial  variation  of  their  structure  is  not  resolved 
adequately  by  the  grid  size  employed  here.  Thus  such  questions  as  shock-wave 
refraction  and  mass  flow  diversion  in  the  vicinity  of  the  sonic  circle  must 
await  further  calculations,  done  with  greater  resolution. 
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